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TITLE: LMMSE-BASED RAKE RECEIVER WITH CHANNEL TAP ASSIGNMENT 

BACKGROUND — FIELD OF APPLICATION 

This invention relates to receiver design in communications systems, specifically to improv- 
ing the system performance in terms of the bit error rate (BER) and frame error rate (FER), and 
to reducing the amount of computations and/or hardware in an effective and reliable manner. 

BACKGROUND — DESCRIPTION OF PRIOR ART 

Conventional Rake Receiver 

A direct-sequence code division multiple access (DS-CDMA) system typically employs a 
rake receiver to recover information sent over a communication medium (channel), e.g., air in 
wireless systems. A conventional rake receiver is based upon the concept of resolvable multi- 
paths. Multipaths are typical of wireless channels where a mobile terminal (cell phone) receives 
multiple copies of the signals from a base station due to reflections from the surroundings. FIG. 
1 illustrates the multiple-path phenomenon in a cellular system. A mobile terminal 20 conmiuni- 
cates with a base station 22 in a cell 24 (base cell). In FIG. 1, there are a signal path 30a directly 
from base station 22, a signal path 30b reflected by a building 26, and a signal 30c reflected by a 
mountain 28. When the delays between the multipaths are large enough, the rake receiver is able 
to distinguish each multipath from others, and the multipaths in such a scenario are referred to as 
"resolvable multipaths". FIG. 2 illustrates a channel impulse response (CER) seen by the rake 
receiver with three well-separated, or resolvable, multipaths. 
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A conventional rake receiver demodulates each multipath separately. In a conventional rake 
receiver, each multipath demodulator is referred to as a "finger". The receiver then combines the 
output of each finger to obtain the demodulated data. 

A conventional rake receiver has following disadvantages: 

1. Intracell interference degrades the output signal-to-noise ratio (SNR) of the rake receiver. 
Herein the SNR refers to both the signal-to-noise ratio and the signal-to-interference- 
plus-noise ratio (SINR) unless stated otherwise. When the rake receiver demodulates 
each multipath, it treats interferences of other multipaths as noise. Such noise is referred 
to as intracell interference. Intracell interference cannot be suppressed by increasing the 
input signal strength, since the strengths of all multipaths increases proportionally to the 
strength of the input signal. As a result, the output SNR of the rake receiver improves 
much more slowly with the increasing input SNR. Eventually the output SNR is limited 
by an SNR ceiling. FIG. 3 illustrates such an output SNR ceiling effect. The intracell 
interference severely limits the BER/FER performance, and the capacity of, e.g., wireless 
cellular systems. 

2. Each finger of the receiver requires a dedicated tracking loop. The rake receiver uses the 
tracking loop to track the delay (timing) and complex amplitude of the multipath assigned 
to a finger. The tracking loops add complexity to the rake receiver. Herein the term 
"complexity" denotes any measure of a receiver with respect to the amount of computa- 
tions, size of the software code, memory size, circuit size, gate count, silicon area, power 
consumption, cost, etc. 

3. It is difficult to track closely spaced multipaths. If multipaths become close to each 
other, the outputs of the tracking loops will contain large errors, as the multipaths are no 
longer well separated as in FIG. 2. The estimation errors in timing and amplitude further 
degrade the BER/FER performance. Sophisticated algorithms can be used to track 
closely spaced multipaths, see, for example, G. Fock, et al, "Channel Tracking for Rake 
Receivers in Closely Spaced Multipath Environments", IEEE Journal on Selected Areas 
in Communications, vol. 19, no. 12, pp. 2420-2431, December 2001, incorporated by 
reference herein. Not only such algorithms significantly increase the complexity of the 
rake receiver, they also require a priori knowledge of information on delays and complex 
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amplitudes for the multipaths to be tracked. Such information may not be available or 
can be hugely erroneous when estimated from closely spaced multipaths, thus degrading 
the performance of the tracking loops. 

4. The quality of the initial multipath estimation is poor when the multipaths are closely 
spaced. As mentioned in previous paragraph, this has adverse effect on tracking algo- 
rithms. This also makes the finger assignment error-prone. 

5. A conventional rake receiver based upon the resolvable multipath concept typically re- 
quires a sample rate that is four times (4x) the chip rate (the symbol rate of the spreading 
sequence) or higher, to reduce the penalty in signal-to-noise ratio (SNR) due to mis- 
alignment of the peak of a multipath and its corresponding sample point. Herein a sam- 
ple rate which is four times the chip rate is referred to as 4x oversampling, and a sample 
rate which is two times the chip rate is referred to as 2x oversampling, and so on. For 
example, 2x oversampling results in about 0.5 dB SNR loss, while 4x oversampling re- 
duces the SNR loss to about 0.1 dB. Higher sampling rate, however, increases the power 
consumption of the analog-to-digital converter (ADC) and front-end pulse-matched filter, 
if it is implemented digitally, which is undesirable especially in a mobile terminal that is 
power-limited. 

Generalized Rake Receiver 

A generalized rake receiver (G-rake) addresses the intracell interference problem in a con- 
ventional rake receiver. For details, see Y.-P. E. Wang, J.-F. Cheng and E. Englund, "The Bene- 
fits of Advanced Receivers for High Speed Data Conmiunications in WCDMA", Proceedings of 
IEEE Vehicular Technology Conference, pp. 132-136, September 2002, incorporated by refer- 
ence herein. In short, G-rake places more fingers than there are multipaths to provide better sup- 
pression of both intracell interference and intercell interference (interference from the base sta- 
tions in other cells). 

A G-rake receiver has following disadvantages: 

1. G-rake receiver uses the multipath decomposition of the CIR, which makes it very sensi- 
tive to channel estimation errors. Let h{t) be the composite CIR, then the multipath 
decomposition of h(t) can be written as 
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Ht) = J^aip(t'-Xi), (1) 

where a/ and X/ are the complex amplitude and delay of the /-th multipath, respectively, 
L is the number of multipaths, and p{t) is the waveform of the CIR when L = 1 . Herein 
the term "composite CIR" refers to the CIR h{t) in its entirety, as opposed to its multi- 
path decomposition on the right hand side of Equation (1). G-rake uses the multipath de- 
composition to compute the noise covariance matrix that is needed in computing the op- 
timum weights for combining the finger outputs. The multipath decomposition requires 
estimation of the complex amplitude and delay of each multipath, and thus is susceptible 
to channel estimation noise, especially when the multipaths are closely spaced. As a re- 
sult, in real-world scenarios G-rake loses much of its promised performance gains that are 
predicted under the perfect knowledge of the channel. Details of the realistic G-rake per- 
formance can be found in, for example, G. Kutz and A. Chass, "On the Performance of a 
Practical Downlink CDMA Generalized Rake Receiver", Proceedings of IEEE Vehicular 
Technology Conference^ pp. 1352-1356, September 2002, incorporated by reference 
herein. 

2. Efficient and effective finger assignment algorithms have yet to be developed for a G- 
rake receiver. The performance of interference suppression depends strongly on the loca- 
tions of additional fingers. The search for optimum locations of additional fingers, how- 
ever, is computationally prohibitive. For example, each search step requires inversion of 
the noise covariance matrix. U.S. patent application 09/845,950 to Y.-P. E. Wang, G. E. 
Bottomley and K. Urabe, discloses an iterative method in lieu of matrix inversion. The 
iterative method, however, is not guaranteed to converge. The solution given by a di- 
verging iteration can severely degrade the receiver performance. The article of Y.-P. E. 
Wang, J.-F. Cheng and E. Englund, cited heretofore, discloses another finger assignment 
strategy that simply puts two additional fingers one chip before and after the known CIR. 
A heuristic search algorithm was proposed in the article of G. Kutz and A. Chass, cited 
heretofore, proposes a heuristic search algorithm ("Kutz and Chass scheme"). These fin- 
ger assignment schemes, while simple, are not optimum, and thus can be ineffective in 
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interference suppression, especially when the multipaths are closely spaced and the chan- 
nel is severely dispersive. 

3. Other disadvantages of a conventional rake receiver, such as the need for finger tracking 
and high sample rate, still exist in a G-rake receiver. 

LMMSE Receiver 

A linear-minimum-mean-square-error (LMMSE) receiver is typically designed based upon 
the composite CIR. A conventional LMMSE receiver has uniformly spaced taps, and the tap co- 
efficients can be computed based on the composite CIR to minimize the mean output error en- 
ergy. An LMMSE receiver is thus robust in the presence of the channel estimation errors. For 
details, see A. Mirbagheri and Y. C. Yoon, "A Linear MMSE Receiver for Multipath Asynchro- 
nous Random-CDMA with Chip Pulse Shaping", IEEE Transactions on Vehicular Technology, 
vol. 31, no. 5, pp. 1072-1086, September 2002. 

A conventional LMMSE receiver has following disadvantages: 

1. The number of uniformly spaced taps is much more than the number of taps (fingers) in a 
rake receiver, thus having more complexity. 

2. To improve the performance, an LMMSE receiver will have fractionally-spaced taps, 
meaning more than one tap per chip. This further increases the complexity, and may be 
numerically unstable when the number of taps is large. 

3. To obtain the optimum tap coefficients requires inversion of the noise covariance matrix. 
The dimension of the covariance matrix is the same as the number of taps, so this be- 
comes impractical for large number of taps. 

4. To avoid matrix inversion, an LMMSE receiver can be made to adapt to changes in CIR. 
The amount of computations in adaptive algorithms is often too huge to be practical. 
Adaptive algorithms also cause loss in performance. 
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SUMMARY OF THE INVENTION 

The present invention relates to receiver design in communications systems, specifically to 
computationally efficient receiver implementations that effectively and reliably improve the 
BER/FER performance of conmiunications systems. 

In accordance with the present invention, an LMMSE-based rake receiver selects channel-tap 
locations based on the composite CIR. The present invention includes two types of channel-tap 
selection schemes: sequential search and heuristic search. A sequential search uses a recursive 
evaluation to efficiently select channel-tap locations, which yields a near-optimum performance 
while providing efficient implementations. A heuristic search has a very simple implementation, 
and in many cases has a performance comparable to that of a sequential search. Both types of 
channel-tap selection schemes provide substantial SNR gain over a conventional rake receiver. 
Use of the composite CIR makes the receiver robust in the presence of the channel estimation 
errors, and eliminates the need for tracking loops and high sample rates. 

OBJECTS AND ADVANTAGES 

Accordingly, several objects and advantages of the present invention are: 

(a) to provide a receiver that has better and more stable BER/FER performance than existing 
rake receivers in the presence of interferences; 

(b) to provide a receiver that maintains its performance in the presence of the channel estima- 
tion noise; 

(c) to provide a receiver that eliminates the need for finger-tracking functions in existing 
rake receivers; 

(d) to provide a receiver that has similar performances at 2x oversampling rate and at 4x 
oversampling rate; 

(e) to provide a receiver that is power-efficient; 

(f) to provide channel tap assignment schemes that select locations for channel taps in an ef- 
ficient and effective manner to achieve optimum or near-optimum performance; 

(g) to provide channel tap assignment schemes that keep the number of the channel taps as 
minimum as possible. 
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Further objects and advantages will become apparent from a consideration of the ensuing de- 
scription and drawings. 

BRIEF DESCMPTION OF THE DRAWINGS 

The present invention is best understood from the following detailed descriptions when read 
in connection with the accompanying drawings. In the drawings, like numbers refer to like ele- 
ments. Included in the drawings are the following figures: 

FIG. 1 illustrates the multipath phenomenon in a wireless cellular system. 

FIG. 2 illustrates a well-separated multipath channel impulse response seen by a rake re- 
ceiver. 

FIG. 3 illustrates the ceiling effect on the output SNR of a conventional rake receiver. 

FIG. 4 illustrates an exemplary LMMSE-based rake receiver according to certain embodi- 
ments of the present invention. 

FIG. 5 illustrates an exemplary non-uniformly spaced LMMSE receiver according to other 
embodiments of the present invention. 

FIG. 6 illustrates an exemplary flowchart of sequential search operations of selecting chan- 
nel-tap locations according to certain embodiments of the present invention. 

FIG. 7 illustrates another exemplary flowchart of sequential search operations of selecting 
channel-tap locations according to other embodiments of the present invention. 

FIG. 8 illustrates an exemplary flowchart of operations of pre-selecting channel taps accord- 
ing to embodiments of the present invention. 

FIG. 9 illustrates a contiguous search region with respect to the span of the CIR according to 
embodiments of the present invention. 

FIG. 10 illustrates a search region consisting of four disjoint sub-regions for a well-separated 
two-multipath channel according to one embodiment of the present invention. 

FIG. 11 illustrates a search region consisting of two disjoint sub-regions for a well-separated 
two-multipath channel according to another embodiment of the present invention. 
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FIG. 12 illustrates an exemplary flowchart of heuristic search operations of selecting chan- 
nel-tap locations according to embodiments of the present invention. 

FIG. 13 illustrates an exemplary flowchart of operations of selecting channel-tap locations in 
a heuristic search scheme according to certain embodiments of the present invention. 

FIG. 14 illustrates an exemplary flowchart of operations of selecting channel-tap locations in 
another heuristic search scheme according to other embodiments of the present invention. 

DETAILED DESCRIPTION 

The present invention relates to the receiver design in conmiunications systems. In a com- 
munications system, the receiver has means of estimating the composite CIR h(t) . For example, 
in a wireless CDMA system, h(t) can be estimated from a pilot channel. Given the channel-tap 

locations ro»^i» - '^M-i' let h = [hQ,hi,,.,,hj^_if = [M^oX'^C^iX- -M^m-i)]^ be the channel 
response vector. The received vector r with respect to the data symbol b can be written as 

r = hZ7-hz, (2) 

where z is the noise vector. In a wireless CDMA cellular system, z can be modeled as the sum of 
the intracell interference and additive white Gaussian noise (AWGN), where the AWGN consists 
of the intercell interference and thermal noise. Consequently, the covariance matrix R of z can 
be written as 

^ ^BC + ^AWGN » (3) 

where Rg^ ^ awgn contributions to R from the intracell interference (the interfer- 

ence from the base cell, the base-cell component) and the AWGN, respectively. 

The (i, j)-ih element of R^c is proportional to the following quantity: 

cov(z,.,z;)bc =^BcZ^(^' +'^^c)^*(0 ■^'^^c)^ (4) 

where E^q is the total received signal energy from the base cell (desired user signal energy + 
other user signal energy), and is the chip duration. Estimation of E^q will be described later. 
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Equation (4) is derived for a wireless CDMA cellular system with orthogonal spreading 
codes. For systems using random, non-orthogonal spreading codes, Equation (4) becomes 

cov(z,.,Zy)Bc =^BcZ'^(^/ + nT^)h\t j nT^) - Eohiti)h\t j) , (5) 

n 

where Eq is the signal energy of the desired user. 

The (/, j>th element of Ry^^c^ proportional to the following quantity: 

C0V(Z,- , AWGN = ^og(h - 0* ) ' 

where Nq is the single-sided power spectrum density of the thermal noise, and g(t) is the auto- 
correlation function of the impulse-response function of the front-end baseband filter. Estima- 
tion of Nq will be described later. 

It should be noted that the covariance expressions in Equations (4), (5) and (6) may differ 
from those used in various implementations by a constant scaling factor. The scaling factor is 
determined by the spreading factor, signal strength, path gains, etc. 

The LMMSE receiver chooses a weight vector w so that the quantity MSE = E[\w^r -b\^] 
is minimized. The superscript H denotes the Hermitian transpose. Such a weight vector is given 
by 

w^ocR^^h (7) 
where a is a scaling factor. The MMSE is proportional to the following quantity: 

, H , > (8) 

where P is a scaling factor. 

The weight vector can be also calculated to optimize other design criteria, such as maximum 
signal-to-interference-plus-noise-ratio (SINR). For linear receivers, the MMSE and maximum- 
SINR criteria give the same results. 

FIG. 4 illustrates an exemplary receiver in accordance with certain embodiments of the pre- 
sent invention. A channel estimator 402 estimates the composite CIR from the received signal. 
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A tap selector 404 selects channel-tap locations based on the composite CIR provided by channel 
estimator 402. Herein the term "channel-tap location" and the term "channel tap" are used inter- 
changeably so long as the meaning is clear in the context. A weight calculator 406 computes the 
weight vector using the output of tap selector 404. A data demodulator 408 recovers the infor- 
mation from the received signal using the channel-tap locations from tap selector 404 and the 
weight vector from weight calculator 406. In a wireless CDMA cellular system, data demodula- 
tor 408 consists of a correlator with the correlation timing determined by tap selector 404, and a 
combiner with the combining coefficients provided by weight calculator 406. Note that a corre- 
lator may consist of one to several physical hardware correlators. 

The receiver in FIG. 4 is referred to as an LMMSE-based rake (L-rake) receiver. Typical 
applications of an L-rake receiver are wireless CDMA systems, just like a conventional rake re- 
ceiver. An L-rake resembles a conventional LMMSE receiver in that (i) the weight vector in an 
L-rake receiver is designed based on the MMSE criterion and (ii) an L-rake receiver requires 
only the composite CIR instead of its multipath decomposition, which makes its performance 
very robust in the presence of the channel estimation errors. 

An L-rake receiver also resembles a conventional rake receiver in that the channel taps 
(called "fingers" in a conventional rake receiver) are assigned according to the channel condi- 
tions, as opposed to the uniformly spaced channel taps in a conventional LMMSE receiver. 
When the channel does exhibit well-separated multipaths, the channel taps that would be selected 
by a conventional rake receiver will likely be selected by an L-rake receiver as well. An L-rake 
receiver, however, abandons the concept of resolvable multipaths altogether, and thus does not 
need finger-tracking loops. 

In accordance with other embodiments of the present invention, a non-uniformly spaced 
LMMSE receiver replaces a data demodulator with an LMMSE filter. Such a receiver is a more 
general version of an L-rake receiver, and can also be used in other applications than wireless 
CDMA systems. FIG. 5 illustrates such a receiver. Channel estimator 402 estimates the com- 
posite CIR from the received signal. Tap selector 404 selects the filter-tap locations based on the 
composite CIR provided by channel estimator 402. Weight calculator 406 computes the filter 
coefficients using the output of tap selector 404. LMMSE filter 420 filters the received signal 
with the filter-tap locations determined by tap selector 404 and filter coefficients determined by 
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weight calculator 406. The LMMSE receiver in FIG. 5 can operate at any fractional oversam- 
pling rate. A fractional oversampling rate refers to a sample rate that can be non-integer multiple 
of the symbol rate. Thus in addition to integral oversampling rates such as 2x or 4x oversam- 
pling, a fractional oversampling rate can be 3/2x oversampling, 4/3x oversampling, etc. 

Channel Tap Assignment - Sequential Search 

Given the total number of channel taps Nj , an optimum channel tap assignment would 
search over all combinations of Nj channel-tap locations and choose the one that optimizes a 
pre-determined design criterion. Such an optimum search, however, is far from practical because 
of the huge number of combinations of the channel-tap locations, and amount of computations 
incurred. In contrast, a sequential search selects channel-tap locations one at a time. In a se- 
quential search, a new channel tap is chosen that optimizes a design criterion given all previously 
chosen channel-tap locations. Thus the optimization in a sequential search is one-dimensional, 
as opposed to a multi-dimensional optimization problem in a full-blown optimum search. 

A sequential search scheme for selecting the channel taps, which provides almost identical 
performance as an optimum search but can be implemented very efficiently, will be now de- 
scribed as follows. The sequential search begins with determining a search region. A search re- 
gion is typically a contiguous region beyond which the energy of the composite CIR is negligi- 
ble. Methods of determining the search region will be described in later. 

The search is conducted over all sample points within the search region at, e.g., 2x or 4x 
oversampling rate. Assuming that there are n channel-tap locations that have been selected: 

hA^'"^^n-i- =[^0'^i'--"^«-il^ =[KtQXhiti),..,,hit^_i)f and R„ bethecorre- 

sponding noise covariance matrix. The objective here is to find the next tap /i„ = h(t^ ) that 

minimizes the quantity in equation (8). From Equation (8), the MMSE is minimized when the 
quantity 

Y„+i =h"+iR;;iih„+i (9) 
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is maximized. Direct evaluation of Equation (9) requires inversion of the covariance matrix for 



each candidate channel tap. Matrix inversion of dimension n involves 0(n ) operations. A 



more efficient approach is to use the following recursive evaluations: 
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and 
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where a J =cov(z„,z„) and r„ =[cov(zo,Z;z),cov(zi,z„), ••,cov(z„_i,z„)]^ is called the new 
covariance vector. Recursive equations (10) and (11) make it explicit the dependence of y„+i 

and R~|i on the previously evaluated functions of previously chosen channel-tap locations and 
on the functions of the next channel-tap location r„ . 

To maximize Y^+i , the sequential search finds the next channel tap = h(t^ ) to maximize 
the quantity 
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in Equation (10) since Yn is independent of . Evaluation of Equation (12) requires 0(n ) 

operations. The recursive evaluation begins with Rq^ =Oq =cov(zo>^o) Yi =|'^o P ^^o 
The recursive evaluation avoids matrix inversion, and yet provides accurate results. 

Using an approximate recursive evaluation may result in a still more efficient sequential 
search scheme. For example, the quantity R^^h„ in Equation (12) is independent of the new 
channel tap and so can be pre-computed for use in the evaluation of each candidate tap. If the 
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term r"R„ t„ in Equation (12) is dropped, the next channel tap can be chosen to maximize the 
quantity 

I' , (13) 

which requires only 0(n) operations. The performance difference between the exact sequential 
search and the approximate sequential search using Equation (13) is negligible in most cases. 

The sequential search stops when the number of selected channel taps reaches the total num- 
ber of taps Nj , which is typically determined by the hardware and other design constraints. 
Alternately, after each new channel tap is selected, the sequential search checks whether the new 
channel tap significantly improves the performance. If the performance is not improved noticea- 
bly, indicating more additional channel taps will be unlikely to generate large performance gains, 
the sequential search stops early, so the receiver uses as few channel taps as possible for satisfac- 
tory performance. The quantity in Equation (12) can be used for measuring the performance im- 
provement. The sequential search can also stop early if Equation (10) reveals that a predeter- 
mined value of the design criterion has been met. 

The sequential search can start with n = 0. Alternately the search can start with n = N^ with 
the first channel taps being pre-selected. This further reduces the amount of computations 
without degrading the performance, if the first A^s taps have strong energy and are well sepa- 
rated. In particular, choosing the first tap with the maximum energy always maximizes, or 
nearly maximizes Yi • 

FIG. 6 illustrates an exemplary flowchart of sequential search operations within tap selector 
404 in FIG. 4 in accordance with certain embodiments of present invention. Operations in FIG. 
6 start with step 502. Step 502 determines search region 5. At step 504, a test compares if the 
number of channel taps to be pre-selected, , is zero. cm be determined, and supplied to 
the receiver, by the control functions of a conmiunications system where the receiver is embed- 
ded. If the test at step 504 determines that is zero, operations proceed to step 508. If, how- 
ever, the test at step 504 determines that is not zero, operations proceed to step 506. Step 
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506 pre-selects channel-tap locations, fo»^i r *>^yVs-i > the span of the CIR. Herein the 
term "the span of the CIR" refers to a contiguous region Tqu^ beyond which the energy of the 
composite CIR is negligible. Note that may be updated at the end of step 506, depending on 
the actual number of pre-selected channel taps by step 506. 

Step 508 sets a tap counter nto N^. Step 510 checks if all sample points in search region S 

have been searched. If step 510 determines that all sample points in search region 5 have been 
searched, operations proceed to step 522. If, however, step 510 determines that not all sample 
points in search region S have been searched, operations proceed to step 512A. Step 512A se- 
lects tap n that maximizes Equation (12). Step 514 updates recursive equations (10) and (1 1) 
according to the result of step 512A. Step 516 increments the tap counter n by 1. At step 518, a 
test compares if the tap counter n has reached Nj . If the test at step 518 determines that the tap 
counter n reaches Nj , operations proceed to step 522. If, however, the test at step 518 deter- 
mines that the tap counter n has not reached Nj , operations proceed to step 520. At step 520, a 
test compares if the difference in performance measurements, Yn+i " Y« » exceeds a threshold 
Dp . If the test at step 520 determines that the performance difference exceeds the threshold 
Dp , operations proceed to step 510. If, however, the test at step 520 determines that the per- 
formance difference does not exceed the threshold, operations proceed to step 522. Step 522 
outputs the number and locations of the selected channel taps. The threshold Dp will have no 
effect if it is set to be less than zero. 

FIG. 7 illustrates another exemplary flowchart of sequential search operations within tap se- 
lector 404 in FIG. 4 in accordance with other embodiments of present invention. Operations in 
FIG. 7 are identical to the operations in FIG. 6, except that step 512B in FIG. 7 replaces step 
512A in FIG. 6. Step 512B selects tap n that maximizes Equation (13). 

The sequential search schemes in FIG. 6 and FIG. 7 can also incorporate additional algo- 
rithmic constraints. For example, the tap selection in step 512A or 512B can be restricted only to 
the locations the minimum distance of which to any previously selected channel-tap location is at 
least D. The reason for such a restriction is that closely spaced channel taps are not effective in 
suppressing the interference, and is highly sensitive to numerical errors. Thus the restriction 
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helps in reducing the size of the sample points to be searched, and in improving numerical stabil- 
ity, especially at high oversampling rates, such as at 4x oversampling. In embodiment of the 
present invention, D is chosen to be an integer multiple of the sample intervals, such as or 

T^I2 at 2x or 4x oversampling. 

Preselection of Channel-Tap Locations 

In some cases it is desirable to pre-select some channel-tap locations. The objective is to pre- 
select a number of strongest channel taps that satisfy certain distance constraint. Herein a strong 
channel tap refers to a channel tap with large energy. 

FIG. 8 illustrates an exemplary flowchart of channel-tap pre-selection operations within step 
506 in FIG. 6 in accordance with embodiments of present invention. Operations in FIG. 8 start 
with step 602. Step 602 sorts the composite CIR samples by their energy levels. Step 604 ini- 
tializes a CIR sample counter m and a tap counter n to zero. Step 606 picks the tap location 

so that h{t^ ) has the w-th largest energy. At step 608, a test compares if the energy of h{t^ ) 
exceeds an energy threshold . Step 608 ensures that pre-selection stops when there is no 
available candidate channel tap with large energy. 

If the test at step 608 determines that the energy of h{t^ ) does not exceed the energy thresh- 
old Ej , operations proceed to step 620. If, however, the test at step 608 determines that the en- 
ergy of h{t^ ) exceeds the energy threshold , operations proceed to step 610. At step 610, a 

test examines if is within the distance D* of any previously selected channel-tap locations. 

Step 610 guarantees that the minimum distance between the pre-selected tap locations is at least 
D\ In embodiments of the present invention, is chosen to be an integer multiple of the sam- 
ple intervals. D' will have no effect on pre-selection of channel-tap locations if it is set to be 
zero or less. 

If the test at step 610 determines that is within the distance D' of a previously selected tap 

location, operations proceed to step 616. Step 616 increments the CIR sample counter m by I. 
At step 618, a test compares if the CIR sample counter m has reached to M, the total number of 
CIR samples. If the test at step 618 determines that the CIR sample counter m has reached to M, 
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operations proceed to step 620. If, however, the test at step 618 determines that the CIR sample 
counter m has not reached to M, operations proceed to step 606. 

If, however, the test at step 610 determines that is not within the distance D' of any of the 
previously selected channel-tap locations, operations proceed to step 612. Step 612 increments 
the tap counter n by 1, admitting as the new tap location. At step 614, a test compares if the 

tap counter n reaches . If the test at step 614 determines that the tap counter n has not 
reached A^s ' operations proceed to step 616. If, however, the test at step 614 determines that the 
tap counter n reaches A^s , operations proceed to step 620. 

Step 620 updates to be the actual number of the pre-selected tap locations. Step 622 out- 
puts and the pre-selected channel-tap locations. 

Note that if A^g = 1 , the receiver will simply choose a channel tap with largest energy. Some 
steps in FIG. 8, e.g., sorting, can be skipped. 

Determination of the Search Region 

In accordance with embodiments of the present invention, a receiver selects a search region 
according to the span of the CIR. Typically a search region is a contiguous region that includes 
the span Tqi^ of the CIR. As defined previously, Tqi^ is a contiguous region beyond which the 
energy of the composite CIR is negligible. FIG. 9 illustrates a search region with respect to 
^ciR according to certain embodiments of the present invention. FIG. 9(A) and FIG. 9(B) 

make clear the definition of T^jr for various CIR shapes. More specifically, T^r starts from 

the location where the energy of the composite CIR first become apparent, and ends at the earli- 
est location after which the energy of the composite CIR is negligible. The energy of the com- 
posite CIR can be strong everywhere within T^jr , as shown in FIG. 9(A). The energy of the 

composite CIR can be also zero or negligible in some part of T^r , as shown in FIG. 9(B). An 

energy threshold can be used to determine whether the energy of the composite CIR at a sample 
point is considered to be negligible. 
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In FIG. 9, a contiguous search region S consists of three sections: The span of the CIR Tqi^ , 
the pre-CIR section Tl and the post-CIR section . In sections Tl and , the composite 
CIR is zero or has very small amplitudes. The channel taps in Tl and help in reducing the 
interferences. The sizes of and generally depends on the shape of the composite CIR. 
Larger sizes of and result in better performance. In practical situations, however, most 
of the performance gain can be obtained with 7^ and being no larger than Tqi^ . 

If the composite CIR exhibits well-separated multipaths, the size of the search region can be 
reduced. FIG. 10 illustrates a search region for a two-multipath CIR according to one embodi- 
ment of the present invention. In FIG. 10, each path has its own search region referred to as 
"path region". Path region 5| includes path 1 of the composite CIR and path region S2 includes 
path 2 of the composite CIR. Region 53 is the mirror image region of path 2 with respect to 
path 1. Region 54 is the mirror image region of path 1 with respect to path 2. The total search 
region 5 is the union of regions 5^ , ^2 , ^3 and ^4 . 

The mirror region of a path with respect to another path (reflecting path) can be determined 
in several ways. A typical choice is to place the mirror region of a path in a location that is 
symmetrical to the strongest channel tap in the reflecting path. 

While FIG. 10 illustrates a search region for a two-multipath CIR, a search region for a 
more-than-two-multipath CIR can be constructed similarly. In this case, the total search region S 
is the union of all path regions and mirror image regions. 

FIG. 11 illustrates a search region for a two-multipath CIR according to another embodiment 
of the present invention. In FIG. 11, the total search region S consists of path regions and 52 
only. 

Channel Tap Assignment - Heuristic Search 

A heuristic search for selecting the channel taps is an alternative to the sequential search de- 
scribed previously. A heuristic search first pre-selects the channel taps in the span of the CIR to 
capture the signal energy. The heuristic search then proceeds to place additional channel taps at 
certain distances of the pre-selected taps to suppress the interference. 
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FIG. 12 illustrates a flowchart of operations of an exemplary heuristic search. Operations in 
FIG. 12 start with step 506. Step 506 pre-selects channel-tap locations, ^c^i' "'^/Vs-i ' 

the span of the CIR. Step 530 uses a heuristic search scheme to select up to channel taps. 
Step 532 outputs the number and locations of the selected channel taps. 

FIG. 13 illustrates an exemplary flow chart of operations in a heuristic search scheme within 
step 530 in FIG. 12 according to certain embodiments of the present invention. Operations in 
FIG. 13 start with step 630. Step 630 initializes a tap counter n, a pre-selected tap counter p and 
a distance counter q to zero. Step 632 makes a Q-entry list of distinct distances » » ' * » ^e-i 

among the pre-selected channel-tap locations. Step 634 considers if tp - is within the dis- 
tance D" of any already-chosen channel-tap location. If step 634 determines that tp - is 

within the distance D" of an already-chosen channel-tap location, operations proceed to step 
640. If, however, step 634 determines that tp - d^ is not within the distance D" of any already- 
chosen channel-tap location, operations proceed to step 636. In embodiments of the present in- 
vention, D" is chosen to be an integer multiple of the sample intervals. 

Step 636 takes tp -dg as a new channel-tap location, and increments the tap counter n by 1. 

At step 638, a test compares if the tap counter n reaches , a predetermined number of al- 
lowed channel taps for the heuristic search scheme. If the test at step 638 determines that the tap 
counter n reaches , the heuristic search is complete. If, however, the test at step 638 deter- 
mines that the tap counter n has not reached N^i , operations proceed to step 640. 

Step 640 considers if tp + d^ is within the distance D" of any already-chosen channel-tap 

location. If step 640 determines that tp + d^ is within the distance D" of an already-chosen 

channel-tap location, operations proceed to step 646. If, however, step 640 determines that 
tp + dg is not within the distance D" of any already-chosen channel-tap location, operations 

proceed to step 642. 

Step 642 takes tp -td^ as a new channel-tap location, and increments the tap counter n by 1. 
At step 644, a test compares if the tap counter n reaches . If the test at step 644 determines 
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that the tap counter n reaches , the heuristic search is complete. If, however, the test at step 
644 determines that the tap counter n has not reached Ny^ , operations proceed to step 646. 

Step 646 increments the distance counter ^ by 1. At step 648, a test compares if the distance 
counter q reaches Q, the number of entries in the distance list. If the test at step 648 determines 
that the distance counter q has not reached Q, operations proceed to step 634. If, however, the 
test at step 648 determines that the distance counter q reaches Q, operations proceed to step 650. 

Step 650 increments the pre-selected tap counter p by 1. At step 652, a test compares if the 
pre-selected tap counter p reaches . If the test at step 652 determines that the pre-selected tap 

counter p reaches , meaning that all pre-selected channel-tap locations have been considered 

by the heuristic search scheme for selecting additional channel-tap locations, the heuristic search 
is complete. If, however, the test at step 652 determines that the pre-selected tap counter p has 
not reached , operations proceed to step 654. Step 654 resets the distance counter q to zero, 
and operations proceed to step 634. 

The channel-tap locations selected by the heuristic search scheme in FIG. 13 have a property 
that the distance between a channel tap selected by the heuristic search scheme in FIG. 13 and a 
channel tap pre-selected by step 506 in FIG. 12 equals to the distance between a channel tap pair 

pre>selected by step 506 in FIG. 12. This is a special case of a more general heuristic search 
scheme by which the distance of a selected channel tap to any other channel tap, selected or pre- 
selected, equals to the distance between a pair of pre-selected channel taps. 

FIG. 14 illustrates an exemplary flow chart of operations in another heuristic search scheme 
within step 530 in FIG. 12 according to other embodiments of the present invention. Operations 
in FIG. 14 start with step 660. Step 660 initializes a tap counter n, a mirror image counter m and 
a pre-selected tap counter p to zero. Step 662 chooses a tentative channel-tap location 2tp - , 

the mirror image of a channel-tap location with respect to a channel-tap location tp . Step 

664 considers if 2tp - is within the distance of any already-chosen channel-tap location. 

If step 664 determines that 2tp - is within the distance D" of an already-chosen channel-tap 

location, operations proceed to step 670. If, however, step 664 determines that It^ - is not 
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within the distance D" of any already-chosen channel-tap location, operations proceed to step 
666. 

Step 666 takes 2tp -t^ as a new channel-tap location, and increments the tap counter n by 

1. At step 668, a test compares if the tap counter n reaches , a predetermined number of al- 
lowed channel taps for the heuristic search scheme. If the test at step 668 determines that the tap 
counter n reaches Nyi , the heuristic search is complete. If, however, the test at step 668 deter- 
mines that the tap counter n has not reached iV^ , operations proceed to step 670. 

Step 670 increments the mirror counter m by L At step 672, a test compares if the mirror 
counter m reaches . If the test at step 672 determines that the mirror counter m has not 

reached , operations proceed to step 662. If, however, the test at step 672 determines that the 
mirror counter m reaches , meaning that all image locations of pre-selected channel-tap loca- 
tions with respect to the pre-selected channel-tap location tp have been considered, operations 
proceed to step 674. 

Step 674 increments the pre-selected tap counter p by 1. At step 676, a test compares if the 
pre-selected tap counter p reaches , If the test at step 676 determines that the pre-selected tap 
counter p reaches A^s ' meaning that all pre-selected channel-tap locations have been considered 
by the heuristic search scheme for selecting additional channel-tap locations, the heuristic search 
is complete. If, however, the test at step 676 determines that the pre-selected tap counter p has 
not reached , operations proceed to step 678. Step 678 resets the mirror counter m to zero, 
and operations proceed to step 662. 

The channel-tap locations selected by the heuristic search scheme in FIG. 14 have a property 
that any channel tap selected by the heuristic search scheme in FIG. 14 is the mirror image of a 
channel tap pre-selected by step 506 in FIG. 12 with respect to another channel tap pre-selected 
by step 506 in FIG. 12. This makes the heuristic search scheme in FIG. 14 a special case of the 
heuristic search scheme in FIG. 13. 

The heuristic search scheme as appeared in FIG. 12 and FIG. 14 is similar in form to the 
Kutz and Chass scheme, in assigning additional channel taps. The main difference is that the 
present invention pre-selects the channel taps based solely upon the energy levels of the compos- 
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ite CIR channel taps, and thus discards the conventional practices of assigning channel taps ac- 
cording to the locations of the multipaths. This brings forth the benefits of not requiring the fin- 
ger tracking, as well as the ability of operating on 2x oversampling. Another difference is that 
the present invention allows the minimum distance D" among the channel taps to be less than the 
one-chip minimum distance specified in Kutz and Chass scheme. 

Modeling Intercell Intetference 

When the intercell interference, i.e., the interference from other surrounding cells, is strong, 
the performance can be improved if the intercell interference is modeled explicitly in the noise 
covariance. In such a case. Equation (3) becomes 

^ - ^BC + ^^OC '^AWGN » (14) 

where Rqc ^ awgn ^he contributions to R from the intercell interference (the interfer- 
ence from other cells, the other-cell component) and the thermal noise (the thermal-noise com- 
ponent), respectively. 

The (/, j>th element of Rqc is proportional to the following quantity: 

p-i 

COV(z,.,Zj)oc = S^0C,pZ^0C,/7(^i +«^c)^OC,p(0 +"^c)' (15) 
p=Q n 

where Eqq ^ and /^ocp (0 are the interference energy and the composite CIR of the ;7-th sur- 
rounding cell, respectively, and P is the number of surrounding cells, /^ocp (^) t)e estimated 
from the pilot signal of the p-th surrounding cell. Estimation of £^oc,p ^^'^ described later. 

It should be noted that in typical wireless networks, the intercell interference is dominated by 
a very few, and often by only one, surrounding cells. The interference from non-dominating 
cells can be modeled as AWGN without causing performance loss. This reduces the number of 
the composite CIRs that need to be estimated and simplifies evaluation of the noise covariance 
matrix. When the intracell interference is dominant, all intercell interference can be modeled as 
AWGN without significant performance loss. 
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Pre-calculation of Noise Covariance Terms 

During a sequential search, the new covariance vector r„ is computed at each candidate tap, 

which involves evaluation of the covariance components defined in Equations (4), (6) and (15). 
For efficient implementations, the terms in those covariance equations can be pre-computed and 
stored so most of the covariance computations can be replaced by simple table accesses. The 
noise covariance can be decomposed into a one-dimensional part, a cyclostationary part, and a 
two-dimensional part. The covariance terms that can be pre-computed are described as follows. 

Rewrite Equation (4) as 



cov(z,-, Zj)^Q = £bc Z + '^T^)h\tj + nT^) 



BC 



V n 



Define 



(16) 



(17) 



and 



(18) 



Functions r^it^Jj) and rQQpit^jj) are cyclostationary in that rf^it|Jj) = r^{ti+T^,tj+T^) and 

rQQ p (ti Jj) = ^oc, p ' ^7 ) • The cyclostationary property allows each of (f,- , tj ) and 

rQQ^pifi.tj) to be stored in several one-dimensional tables. For example, at 2x oversampling, 

ri^(tiJj) can be stored in two one-dimensional tables, and at 4x oversampling, r^it^^tj) can be 

stored in four one-dimensional tables, and so on. Equation (6) requires a single one-dimensional 
table to store one-dimensional function g{t) . The last term of Equation (16) remains two- 
dimensional, but it does not need to be pre-computed as it involves only one simple multiplica- 
tion. 
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Evaluation of the noise covariance requires the knowledge of the composite CIR h{t) . In re- 
gions where the amplitude of h(t) is small, h(t) can be approximated by zero. This reduces the 
length of the non-zero region of h(t) , so the amount of computations is further reduced. 

Energy Level Estimation 

Estimation of the signal and noise energy levels ^bC' ^oc,p ^^^^ t)e described as 

follows. Let y{t) be the received signal consisting of signals from the base cell and other cells, 
and the thermal noise. The autocorrelation function (t,x) is given by 

ryy(t,x) = E[y{t-hT)y*(t)] 

= EbcX^(^ + 1^ - )h\t + T - nr^ ) + . (19) 

n 

Z^OCp Z^OCp + + )hoc.p + '^T; ) + Nog{T) 
p=0 n 

^yyOy"^) Can be computed directly from y(t) for a given timing pair (r,T) . Evaluating 
^yy (^^^) on a set of such timing pairs, Equation (19) gives a set of linear equations from which 
^BC' ^OC,p ^0 can be solved for. The dimension of the linear equations should be at least 

P + 2 . Higher dimension (over-determined system) reduces the estimation error and improves 
the numerical stability. 

The autocorrelation in Equation (16) is cyclostationary, and can be approximated by its time- 
average within one chip interval to produce the averaged autocorrelation that is one-dimensional. 

Multiple Base Cells 

When a mobile terminal receives from multiple base stations, such as during a soft handoff, 
the interference and noise in the signals from one base cell can be modeled to be uncorrelated 
with the interferences and noises in the signals from other base cells because each base cell uses 
a different scrambling sequence. The receiver assigns channel taps using a sequential search or a 
heuristic search, and computes the optimum weight vector for each individual base cell as de- 
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scribed before. The outputs of the receiver for all base cells are then combined to minimize the 
overall output MSE or to maximize the overall output SNR. 

Multiple Receive Antennas 

When a receiver is equipped with multiple antennas receiving from a single cell, the interfer- 
ence and noise in the signals from one antenna are generally correlated with the interferences and 
noises in the signals from the other antennas. For the purpose of illustration, consider a receiver 

with two antennas. In this case, the two received signals will have two time-axes t^^^ and t^^^ , 

The channel-tap locations consist of the channel-tap locations in t^^^ axis and the channel-tap 

locations in t^^^ . The noise covariance between two channel-tap locations in the same time-axis 
remains in the same form as Equations (4), (6) and (15) with respect to the time-axis of interest. 
The covariance between two tap locations in different time-axes takes a slightly different form. 
For example, the base-cell component of the covariance can be written as 

coyi4\zf\c=E^ci:h^\4'^^nT,)h*^'\tf^^nT,), (20) 

where h^^^{t^^^) and h^^\t^^^) are the composite CIRs of the two received channels. The other- 
cell component of the covariance can be modified similarly. The AWGN component of the co- 
variance is zero when two channel-tap locations are in different time-axes, assuming that the cir- 
cuit noises in different antenna receive paths are uncorrected. 

The search process now includes assigning channel taps in both time-axis. Alternately the 
search may take a sub-optimal approach, assigning channel taps and calculating the weight vec- 
tor in one time-axis at a time. The results from the two time-axes are then combined to generate 
the final output. Receivers with more than two receive antennas can be implemented in a similar 
fashion. 

Conclusion, Ramifications, and Scope 

Accordingly, it will be seen that the L-rake receiver of this invention combines the advan- 
tages of both LMMSE receivers and conventional rake receivers, without their disadvantages. 
Use of the composite CIR makes the receiver robust to channel estimation errors, so the receiver 
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performance is largely unaffected by imperfect channel estimations. The channel tap assignment 
schemes eliminate the need for the knowledge of multipath locations, thus making the finger- 
tracking function in a conventional rake receiver unnecessary in the present invention. The pre- 
sent invention also eliminates the problem of initial finger assignment in the presence of closely 
spaced fingers. 

A rake receiver typically requires much fewer taps than an LMMSE receiver does, since an 
LMMSE receiver typically uses equally spaced channel taps regardless of the channel condi- 
tions. Fewer taps reduces the receiver complexity and improves the numerical stability. An L- 
rake receiver retains this advantage of the rake receiver. The sequential search in an L-rake re- 
ceiver places channel taps on optimum or close-to-optimum locations, so that the SNR perform- 
ance is kept at optimum or near optimum level that is far better than the performance of a con- 
ventional rake receiver. Recursive evaluations of the design criterion and inverse of the noise 
covariance matrix make the sequential search highly efficient. Heuristic search provides another 
efficient approach to determining channel-tap locations. 

An L-rake receiver works equally well at 2x and 4x oversampling rates since it does not de- 
pend on the knowledge of the multipath positions. Given that a conventional rake or a G-rake 
receiver typically works at 4x or higher oversampling, an L-rake receiver operating at 2x over- 
sampling offers significant reduction in power consumption. 

An L-rake receiver using a sequential search is able to maintain the fewest possible channel 
taps, since the channel tap assignment process can terminate early if the receiver finds a new 
channel tap does not yield significant performance improvement, or if the performance criterion 
has been met with fewer channel taps. Thus an L-rake receiver is able to keep the power con- 
sumption at minimum. 

The advantages of the L-rake receiver allow it to be implemented efficiently in terms of 
hardware size, silicon area, amount of computations and power consumption. A more general 
form of an L-rake receiver, a non-uniformly spaced LMMSE receiver, extends the advantages of 
the L-rake to all applications where an LMMSE receiver is desired. 

The present invention is disclosed in detail with wireless CDMA systems as illustrative ex- 
amples. The principles of the present invention, however, also apply to other applications where 
an LMMSE receiver is employed. 
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The present invention can be embodied in the form of methods and apparatuses for practicing 
those methods. For example, the present invention can be embodied as circuits or digital hard- 
ware. The present invention can also be embodied in the form of program code stored in tangi- 
ble media, such as floppy diskettes, hard drives, CD-ROMs, CD-RAVs, DVDs, memories, or any 
other machine-loadable storage medium, wherein, when the program code is loaded into and 
executed by a machine, such as a general-purpose computer, digital signal processor, microproc- 
essor, microcontroller, or any other processing circuit, the machine becomes an apparatus for 
practicing the invention. 

It is to be understood that the embodiments and variations shown and described herein are 
merely illustrative of the principles of this invention and that various modifications may be im- 
plemented by those skilled in the art without departing from the scope and spirit of the invention 
as expressed in the following claims. 
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